Theory of thermopower in 2D graphene 
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Motivated by recent experiments we calculate the thermopower of graphene incorporating 

the energy dependence of various transport scattering times. We find that scattering by screened 
charged impurities gives a reasonable explanation for the measured thermopower. The calculated 
thermopower behaves as 1/y/ri at high densities, but saturates at low densities. We also find that 
the thermopower scales with normalized temperature T/Tf and does not depend on the impurity 
densities, but strongly depends on the fine structure constant r 3 and the location of the impurities. 
We discuss the deviation from the Mott formula in graphene thermopower, and use an effective 
medium theory to calculate thermopower at low carrier density regimes where electron-hole puddles 
dominate. 
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Thermopower has been used as a powerful tool to 
probe transport mechanisms in metals and semiconduc- 
tors. Often the measurement of resistivity (or conduc- 
tivity) is inadequate in distinguishing among different 
scattering mechanisms and the thermopower can then 
be used as a sensitive probe of transport properties since 
it provides complementary information to resistivity. In 
this Letter we develop a theory for the thermopower of 
graphene with a goal toward elucidating the comparative 
importance of various scattering mechanisms in graphene 
transport properties. 

Recently, the thermoelectric properties of graphene 
have attracted experimental attention [l|, |2|, y, li| • Exper- 
imentally [J, 0, Q the expected change of sign in the ther- 
mopower is found across the charge neutral point as the 
majority carriers change from electrons to holes. Away 
from the charge neutral region the density dependence 
of thermopower behaves as 1/y/n, and exhibits a linear 
temperature dependence in agreement with the semiclas- 
sical Mott formula [5(. As the temperature increases, a 
deviation from Mott formula is reported [l], Q. Existing 
theoretical works on graphene thermopower either use 
an impurity band model or consider the low tempera- 
ture Mott limit In this paper we present a calcula- 
tion of the thermopower of graphene taking into account 
the energy dependent scattering time for various scatter- 
ing mechanisms. Understanding thermopower requires a 
thorough knowledge of the details of energy dependence 
of transport scattering times [3, H, @] . In metals the Mott 
formula is widely used because the Fermi temperature is 
very high (i.e. T <C Tp) and the scattering time is es- 
sentially energy independent leading to a simple linear-in 
temperature form for thermopower which is proportional 
to the energy derivative of the conductivity evaluated at 
the Fermi energy. Mott formula, derived mathematically 
through the Sommerfeld expansion, is only valid at very 
low temperatures, T/T F -C 1. 

We show that scattering by random charged impurity 
centers, which is the main scattering mechanism limit- 



ing graphene conductivity, (l0( also dominates its ther- 
mopower. We show that the temperature dependent 
screening effects [ll[ must be included in the theory 
to get quantitative agreement with existing experimen- 
tal data. We find the effects of short-range scattering 
and phonons to be negligible in experimental tempera- 
ture range (T < 300 A"), allowing us to ignore phonon 
drag contribution. We also find that the calculated 
thermopower scales with T/Tp and manifests no impu- 
rity density (rii) dependence, but depends strongly on 
the impurity location and the dielectric constant of the 
substrate (or equivalently the fine structure constant of 
graphene). The experimentally observed asymmetry be- 
tween electron and hole thermopower is explained by the 
asymmetry in the charged impurity configuration in the 
presence of the gate voltage. We explain the experi- 
mentally observed sign change near charge neutral point 
(Dirac point) with a simple two component model, which 
we explicitly verify using an effective medium theory cal- 
culation ta king into account the inhomogeneous puddle 
formation. [10( 

The ratio of the measured voltage to the temperature 
gradient applied across the sample is known as the See- 
beck coefficient (or the thermopower) and is given by 
Q = VVyVT, where W is the potential difference and 
VT the temperature difference between two points of the 
sample 7] . In linear response approximation for the elec- 
trical current density, j, and thermal current density, Jq, 
we have: j = Z, U E+Z 12 (-VT), j Q = L 21 E + L 22 (-VT), 

where L l i is defined in terms of the I^ a \ i.e., L 11 = 
7,12 Lr(l) £21 _ _1 r(l) 

j( a ) i s given by 



and L 22 = ^y/ (2) . Here, 
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where is the equilibrium Fermi distribution function, 
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r the relaxation time, and y, the chemical potential, g = 
g s g v is the total degeneracy (g s = 2, g v = 2 being the 
spin and valley degeneracies, respectively), and a(e) is 
the energy dependent conductivity of graphene given by 
a(e) = e 2 VpD(e)T(e) /2, where vp is the Fermi velocity 
and D(e) — g\e\/ (2irh 2 v F ) the density of states. From 
the definition of the thcrmopower we have Q — ju, and, 
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a = L l 

Before we calculate the details of the thermopower for 
different scattering mechanisms we first consider the low 
temperature and high temperature behavior of Q(T). At 
low temperatures (T <^Tp, where Tp = Ep/k B ) we can 
express 1^ as 
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where (3 = l/k B T. Thus we have the well-known Mott 
formula Q of thermopower at low temperatures, i.e. 
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If the energy dependence of the relaxation time is unim- 
portant the sign of the thermopower is determined by 
whether the carriers are electrons or holes. Assuming 
the energy dependent scattering time to be r oc e m we 
have the thermopower at low temperatures 
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We note that in general the exponent 'm' has weak tem- 
perature and density dependence since t behaves only as 
an effective power law in energy. Eq. ([¥]) indicates that 
the thermopower can change sign if m < — 1. At high 
temperatures (T ^> Tp) we can express the 1^ with an 
energy dependent scattering time r = tqs" 1 as 
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where and £(z) are the gamma function and Ric- 
mann's zeta function, respectively, from which we find: 
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FIG. 1: Hole thermopower for different energy dependent 
scattering times, r oc e m . For m < — 1 the low tempera- 
ture thermopower becomes negative. Dashed lines show the 
Mott formula for the corresponding scattering times. Electron 
results are the same with an overall negative sign. 



white noise short-range disorder and acoustic phonon [12j 
scattering we can express the scattering time as r(e) = 
Ti/e, and we have 
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(4) Then the thermopower becomes 
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The thermopower contributions from both neutral scat- 
terers and acoustic phonons are exponentially suppressed 
in the low temperature limit, and can be ignored. For un- 
screened charged impurities we have very simple energy 
dependent scattering time [ll[ r(e) = tq6. Then we have 
the following integrals for t <C 1 
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leading to the thermopower 



At high temperatures Q in graphene approaches a limit- 
ing value. In Fig.Q]we show the calculated graphene thcr- 
mopower for different scattering exponents m (t cx e m ) 
as a function of T/Tp. As shown in Fig. [T] the dashed 
lines representing Mott formula agree with the full cal- 
culations for T < 0.2Tp. In addition the calculated ther- 
mopower scales as a function of the normalized temper- 
ature (T/Tp). 

Now we calculate the thermopower in the presence of 
various physical scattering mechanisms. For both neutral 
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The linear relation of thermopower with temperature (or 
Mott formula) holds to relatively high temperatures in 
graphene for unscreened Coulomb impurities. 

As has been demonstrated theoretically and experi- 
mentally the dominant transport mechanism in graphene 
is the screened Coulomb scattering from charged impu- 
rities. The result of Eq. (TH?)) for unscreened Coulomb 
scattering is much higher than the thermopower observed 



3 



in experiments [l|, 0, E| and cannot explain the behavior 
of Q close to the Dirac point. An accurate quantita- 
tive agreement between theory and experiment can only 
be achieved by taking into account the screening of the 
charged impurities and the strong spatial inhomogene- 
ity that these impurities induce close to the Dirac point. 
For the screened charged impurity scattering, the energy 
dependent scattering time r(efc) is given [l0( 
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where 6 is the scattering angle, v%(q) = 
2ire 2 exp(— qd)/(nq) is the Fourier transform of the 
2D Coulomb potential in an effective background lattice 
dielectric constant k and d is the location of charged 
impurity measured from graphene surface. In Eq. (|13[) . 
e(q) = e(q,T) is the 2D finite temperature static 
RPA dielectric (screening) function appropriate for 
graphene [IlT|. given by e(q,T) = 1 + v c (q)IL(q,T), where 
II(<7, T) is the graphene irreducible finite-temperature 
polarizability function and v c (q) is the Coulomb inter- 
action. There is an important direct T dependence of 
thermopower, not captured in the Mott formula, arising 



from the temperature dependent screening. ll| The 
temperature dependent conductivity due to screening 
effects decreases quadratically at low temperatures 
This mechanism produces a thermopower quadratic 
in temperature rather than linear as in the simple 
Mott formula. Thus we predict a nonlinear quadratic 
temperature correction in the graphene thermopower 
compared with the linear Mott formula. 

In Fig. [2](a) we show the calculated thermopower of 
holes in graphene due to different scattering mecha- 
nisms. The thermopower due to screened charged impu- 
rity is about half of that due to the unscreened charged 
impurity, and increases in a concave manner due to 
temperature and energy dependent screening. On the 
other hand, the thermopower due to neutral scatter- 
ers is exponentially suppressed in the low temperature 
regime. Fig. |2£b) shows the calculated thermopower for 
the screened charged impurity scattering as a function 
of density for different temperatures. As we expect the 
density dependence shows 1/y/n behavior at high den- 
sities. But this power law behavior breaks down and 
saturates at low densities. The saturation value (Q s ) 
does not depend on temperature. The theoretical Q s is 
just a function of interaction (fine structure) parameter 
r s = e 2 /kHvf and the location of impurities d. In Fig. [5] 
(c) and (d) we show the r s and d dependence of the ther- 
mopower. In general the thermopower increases when the 
substrate dielectric constant (k) increases or the charged 
impurities move away. Thus, we predict that the ther- 
mopower of suspended graphene will decrease compared 
with the thermopower of graphene on a substrate for the 
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FIG. 2: (a) Q as a function of temperature for different scat- 
tering mechanisms. Solid, dashed, and dot-dashed line rep- 
resent unscreened Coulomb, screened Coulomb, and neutral 
scatterers, respectively, (b) Q for the screened charged impu- 
rity scattering as a function of density for different tempera- 
tures, (c) Q for different r s — 0.1, 0.85, 2.2 with d — and 
(d) for different d = 0, 5, 10 A with r s = 0.85. In (e) and (f), 
Q in the presence of screened charged impurities for different 
densities is shown. In (e) we use parameters corresponding to 
graphene on S1O2 with mobility fi — 10 4 cm 2 /Vs, and in (f) 
to suspended graphene with mobility fi = 2 x 10 5 cm 2 /Vs. 



same densities because of the reduction of the dielectric 
constant. In Fig [2] (e) and (f) we show the calculated 
thermopower for systems with different mobilities due 
to screened charged impurity scattering for graphene on 
Si02 substrate and for suspended graphene. 

In recent experiments 0, 0, Q it has been observed 
that close to the Dirac point Q does not follow the 
l/\/n scaling predicted by the Mott formula. The rea- 
son for this deviation is that close to the Dirac point 
the quenched disorder induces strong density fluctuations 
that break up the density landscape in electron-hole pud- 
dles 13, H, 14 , EH . To account for the main features of 
the thermopower close to the Dirac point we can use a 
simple two component model near to the charge neu- 
tral regime in which the electron density, n e and the 
hole density nh depend on the doping n according to 
the phenomenological equations: n e = (n rms + n)/2, 
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n)/2 for \n\ < n r 



where n rms is the 



root mean square of the density fluctuations. We assume 
n rms ~ rii fiol . [3| • In the two component model the ther- 
mopower becomes Q = (L 12 + L]^)/(L]} + L^ 1 ) so that, 
for unscreened charged impurities, we find: 
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Thus, if there is an equal number of electrons and holes 
the thermopower goes to zero, and the overall sign is 
decided by the majority carriers. We find that the ther- 
mopower due to the screened charge impurities shows 
the same behavior as shown by the dashed line in Fig. [3] 
In order to explicitly verify the electron-hole puddle pic- 
ture near the charge neutral point we use the effective 
medium theory, EMT for graphene 15| . The density pro- 
file close to the Dirac point is quantitatively described by 
the Thomas-Fermi-Dirac (TFD) theory [lj. Using the 
TFD results the transport properties of graphene close 
to the Dirac point can be accurately calculated using the 
EMT. Denoting by angle brackets disordered averaged 
quantities, for the diagonal transport coefficients from 
the EMT in 2D we have that the effective medium coef- 
ficients, L" s , are implicitly given by the equation [la ]: 



L u (r) 
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Adapting to 2D thermopower the results presented in 
[rH ]. the effective medium off-diagonal coefficient L 12 is 
given by: 
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As shown in [15j for the local values of L l i we can use 
L u (n(r)) and then using Eqs. (fl"5j) . (fT!))) and the proba- 
bility distribution given by the TFD theory we can cal- 
culate the EMT transport coefficients, and in particular 
the effective medium thermopower Q e g = ^Is/^ls f° r 
graphene at, and away from, the Dirac point. The re- 
sults for Q as a function of n at T = 300 K for r s = 0.8, 
d = 1 nm and 7ij = 10 12 cm -2 are shown by the solid 
lines in Fig. 02 in blue (red) are the results obtained with 
(without) the effect of the exchange term on the density 
distribution, we note that the two-component model is 
excellent in describing the main features of the realistic 
EMT Q(n) close to the Dirac point. 

In conclusion we have developed a complete theory for 
the diffusive thermopower of 2D graphene. Quantitative 
agreement between our theory and existing graphene ex- 
perimental thermopower data is a strong indication that 
the dominant carrier scattering mechanism operational in 
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FIG. 3: Q due to the screened charge impurities as a func- 
tion of density close to the Dirac point for T = 300 K, 
m = 10 12 cm -2 , r s — 0.8 and d = 1 nm, obtained using 
the two-component model and the EMT with and without 
exchange energy. 



2D graphene monolayers is screened Coulomb scattering 
by random charged impurities located in the graphene 
environment. At high densities the Mott formula applies 
well to the measured thermopower because of the high 
Fermi temperature, but it fails in low density limit. We 
explain the sign change of the thermopower in the low 
density regime by using both a two component model 
and a realistic effective medium theory, that correctly 
describes transport in the presence of the strong carrier 
density inhomogeneities, that characterize the graphene 
density landscape close to the Dirac point. We make a 
number of specific predictions for graphene thermopower 
(e.g. nonlinearity in temperature, existence of a satura- 
tion thermopower at low densities, nontrivial dependence 
on the background dielectric constant and on the impu- 
rity location), which should be tested experimentally in 
order to conclusively settle the issue of dominant carrier 
scattering mechanism in graphene. 
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